Europhysics Letters 



PREPRINT 



Nonextensive diffusion as nonlinear response 

James F. Lutsko 1 and Jean Pierre Boon 1 

1 Center for Nonlinear Phenomena and Complex Systems Universite Libre de Bruxelles, 
1050 - Bruxelles, Belgium (E-mail: jlutsko@ulb.ac.be; j pboon@ulb.ac.be ) 



PACS. 02.90.+m - Statistical Physics. 
PACS. 05.60.-k - Transport processes. 

PACS. 05.10.Gg - Stochastic analysis methods (Fokker-Planck equation). 



Abstract. - The porous media equation has been proposed as a phenomenological "non- 
extensive" generalization of classical diffusion. Here, we show that a very similar equation 
can be derived, in a systematic manner, for a classical fluid by assuming nonlinear response, 
i.e. that the diffusive flux depends on gradients of a power of the concentration. The present 
equation distinguishes from the porous media equation in that it describes generalized classical 
diffusion, i.e. with r/yDt scaling, but with a generalized Einstein relation, and with power- law 
probability distributions typical of nonextensive statistical mechanics. 



One of the characteristic features of the nonextensive thermodynamics introduced by Tsal- 
lis is the appearance of non-exponential distribution functions with power-law tails [1,2]. There 
has been considerable interest recently in the question of how such non-exponential distribu- 
tions might arise from first-principle considerations [3-6]. In this paper, we show that it is 
possible to obtain, in a generic manner, power-law distributions for the diffusion of a tracer 
particle in a liquid by assuming a simple generalization of the usual linear response arguments. 

In classical transport theory, the probability to find the diffuser at point r at time t given 
that it starts at point r o at time obeys the advection-diffusion equation 

—P (7*, t; 7*0,0) = _r . t?(7,t) P (7, t; + ^=> ■ £ • (7*, 7* , 0) , (1) 

ot or a r a r 

where ~u(l^,t) is the drift and D the diffusion tensor, see e.g. [7]. Note that in general, 
we could, instead of probabilities, speak of the concentration of a diffusing species or the 
density of a single-component system undergoing self-diffusion just as well. There are currently 
two widely explored generalizations of this classical diffusion equation. The first involves a 
fractional time derivative in the diffusion term [8]. This fractional diffusion equation can 
be derived, e.g., by considering the behavior of a diffuser on a lattice with the waiting time 
between hops being governed by a power-law distribution (e.g., a Pareto distribution), a type 
of Levy flight [8] . A second generalization of the classical diffusion equation sometimes used 
to model anomalous diffusion is the so-called "porous media equation" 

-P (7, t; 7*0, 0) = • H (7\ f) P (7*, t; 7*0, 0) + ^ ■ D ■ -^P^ (7*, t; 7* , 0) . (2) 
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As its name indicates, this equation was established (empirically) for particular physical sys- 
tems [9], and, recently, its formal structure, dtf — A(/ Al ), has been analyzed in the mathemat- 
ical literature [10]; it has also been proposed as a phenomenological " nonextensive" generaliza- 
tion of classical diffusion [11, 12]. In one dimension and in the absence of drift, u ( r ,t) = 0, 
this equation admits of power-law solutions of the form 

where the constant A is fixed by normalization. This function can be written more suggestively 
as 

P(M;0,0) = ^-^e 2 _^-^^^) , (4) 

where the so-called q-exponential is defined by 

e q {x) = (1 + (l-q)x)^ , (5) 

and has the property that lim 9 _^i e q (x) — e x . 

Notice that Eq.© can be viewed as a classical diffusion equation with the effective diffu- 
sion coefficient /j, D P^ 1-1 ( r , t; ~r*o, 0). Since concentration-dependent and density-dependent 
diffusion coefficients are commonly used in the scientific literature, this does not differ radi- 
cally from classical diffusion: what is unusual is that the effective diffusion coefficient vanishes 
when the probability (or concentration) vanishes. 

We consider here a collection of atoms moving under Newtonian dynamics and concentrate 
on describing the probability of finding a tagged particle, say the 0-th particle, at position V 
at time t given it begins at position ~r*o at time t = 0, 

P (?, t; y , 0) = (S (t) -?)5tf (0) - ? )) Q , (6) 

where the position of the tagged particle is q (t) and the notation (...) indicates an ensem- 
ble average over some distribution of initial conditions for the other, non-tagged, particles. 
Differentiation of Eq. JJJJ with respect to time leads to the exact equation 

J^P (V, i; TV 0) = ■ it (?, t) P (V, t; ~ro,0) + t^=> • ~J , (7) 
with the current 

J (?, t- ? , 0) = (it' (t) S (g> (t) (0) - y )) , (8) 

where it ( r , t) is the average macroscopic velocity, it (t) is the instantaneous velocity of the 
tagged particle and its peculiar velocity is it' (t) = v (t) — u (r ,t). The idea of linear re- 
sponse theory is that, aside from convection and in the absence of driving forces or boundary 
conditions, the distribution should relax toward a uniform distribution. Hence, one expects 
the current to be proportional, in some sense, to gradients of the distribution and indeed the 
simplest assumption that J oc VP gives the usual phenomenological Fick's law. Higher order 
terms are expected to involve higher order gradients. Such an expansion is most easily devel- 
oped in Fourier space by noting that the Fourier transform of the current and the distribution 
are 

J , t; r> , o) = (it' (t) exp (ik ■ t (t)) 6 (t (0) - 7> )) o , (9) 
p(M;^o,0) - (exp (ik •-?(*)) 5(t(0)-^o)) Q , (10) 
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and expanding their ratio to get a series of the form 



J /P = ik ■ D (t) 



(11) 



thus giving the desired gradient expansion [13]. In 
given by an Einstein relation, is time dependent on 
the time between atomic collisions, but in the long- 
it tends toward a constant value. 

Comparison of Eqs.© and Q suggests that the 
ized diffusion equation is that the current is no lonj 
the distribution function, but rather by gradients in 
order to generalize the linear-response argument, we 
of the /i-th power of the distribution to be 



general, the diffusion coefficient, which is 
microscopic time scales characterized by 
time limit appropriate to hydrodynamics 

key assumption in the use of the general- 
;er most simply described by gradients in 
the distribution raised to some power. In 
proceed by defining the Fourier transform 



Pu 



and expand in the wavevector to obtain 



(~M;y o ,o) = / d V exp fit ■ ~r) P M ("r,£;"r o ,0) 



J /P„=ik ■ D, 



(12) 



(13) 



which gives 



dV p* 0,0) 



(14) 



where D i (t) is the classical diffusion coefficient which is the same as that occurring in Ea. (|ll|l . 
By expressing ~q (t) in terms of the time-integral of the velocity, the diffusion coefficient 
becomes a function of the velocity auto-correlation function as expected 



Da(t) 



(!?(t')l?'(t)Stf(0)~ro)} o dt' 



dr P» Cr,t;"? ,0) 



(15) 



Then, keeping only the lowest order in (|13fl . we obtain the generalized advection-diffusion 
equation 



0_ 

di 



P ( r\ t; 7*0, 0) = = • u ( r\ t) P (V, t- y , 0) + — - • ^ (t) • -=>P» (V, t; y , 0) , (16) 
or Or Or 



with a generalized Einstein relation, Ea. p5|l . relating the diffusion coefficient to the underlying 
microscopic dynamics and corresponding distribution function. Note that even on hydrody- 
namic time scales, for which D i (t) can be replaced by D \ = lim^oo D i (t), D M (t) remains 
time-dependent due to its dependence on P M ( r , t; "to, 0). 

Scaling solutions of Ea. (|16[l in the absence of drift (u = 0) are easily examined. For 
simplicity, consider the approximation, expected to be valid on hydrodynamic time scales, 
D i (t) ~ lim^oo D i (i) = D ±, and suppose that in (^-dimensions the distribution assumes 
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the form P (T*, t; TV 0) — t da F(r/t a ) for some exponent a (note that the prefactor is dictated 
by normalization). Then, one has 



d~t P» (~r , t; ~? ,Q) = J dT' t- d > m F> 1 {r/t a ) 

= S d r d ^-^ a / F"(x)x d - 1 dx, (17) 



with x = t a \r \ and where S4 is the area of the d— dimensional hypersphere. The generalized 
diffusion coefficient is then _ 

*D „ (t) = ^- 1 )« < D ' (18) 



where 

-1 



5 d / F^x^dx 



(19) 



is a time- independent coefficient. Substitution into Ea. (|16f) gives, after simplification, 



-a (dF{x) + 1? • ^ (a0 ) = ^ • ^ • ^ (x) , (20) 

which, for scaling consistency, is only valid for a = \ . The generalized diffusion equation given 
here therefore describes normal diffusion when there is no drift. However, the distributions 
are not Gaussian. For example, in one dimension the scaling equation becomes 

which is again solved by the canonical g-exponential form as 



f|,| = Jr '""-l-iMM)- 



or, with fi — 2 — q, 



Aq-1 r 2 , 

P (r> t; 0,0) = A r ^ ( 1 - (1 - q) - - - - _ ) 



The mean-squared displacement is then obtained straightforwardly from l|23|) : with 
r/ yJD '^t, one has 

/>oo 

(r 2 (t)) = / drr 2 P{r,t; 0,0) 



r 



r°° { r *2\ — 

-i>; ( jf ^^(i-a-,,^-^-) . ,24, 

Thus, in the absence of drift, Ea. (|16|) describes normal diffusion in the sense that the mean- 
squared displacement grows linearly with time, but it generalizes classical diffusion as the 
distribution (j23(l takes on the form typical of nonextensive systems. 

Although Ea. 1)16(1 is formally the same as Eq.Q, there is an additional self-consistency 
implied by the Einstein relation (|15[1 which raises the question as to whether or not these 
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descriptions can be regarded in any sense as being the same. In some cases they are, be- 
cause Ea. (|16|) is of first order in the time so that any time-dependence generated by the 
self-consistency can be hidden by a non-linear change in the time variable. For example, in 
the absence of drift, let Q ( r , s; r o,0) be a solution to Eq.J2J), written with the variable s 
instead of t, and having initial condition P (~r, 0; 7*0, 0). Then it is easy to see that the solu- 
tion of Ea. (|16[l is P (7, t; "r^o, 0) = Q (7*, s (t) ; 7*0, 0) where the time variables are related by 
the implicit equation 

t(s)= [ ( f dV Q^(r, a; V ,0)) da. (25) 



To illustrate, consider again the case of one dimensional diffusion with no drift term. The 
solution of the porous media equation J5J is given by l@J which, with the replacement t — * s, 
is what is called here Q. Substitution into Ea. l25f) gives 

t(s)= (^J d ~r Q M (rV; 7*0,0)) da = Cs^ , (26) 

for some constant C , thus recovering Eq. l|22|l . 

It is interesting that Ea. l|16|) is very similar to the equation recently derived by Abe and 
Thurner [6] where the derivation begins with the classical picture based on a random walker 
which jumps from position r + A at time t to position r at time t + t with a jump probability 
<j) (A) = 4> (—A), so that the probability that the walker be at r at time t + t is given by 

V{r,t + T)= Jv(r + A,t)<t> (A) dA . (27) 

They consider the more general case in which the probability V (r, t) occuring in the integral 
is replaced by the so-called escort probability 

^("••"=tIto' (28) 

with the result that the diffusion equation becomes 

(r, t) = D v (t) (V, t) + i {T (r, t) - V (V, t; )) , (29) 

with D v [i) defined as in Ea. (|15f) . Equation l|29|) is somewhat difficult to interpret as its 
derivation involves taking the limit r — > in which case the second term on the right is 
problematic. Apart from this problem, the result is strikingly similar to our nonlinear-response 
result including the definition of the time-dependent diffusion coefficient (see Eq.(15) in [6]). 
Comparison of the results also indicates that "nonextensive" expressions can arise without 
the explicit introduction of escort probabilities. 

In this paper, we have developed a method whereby a statistical mechanical derivation 
yields typically nonextensive expressions by generalizing a key ansatz, i.e. a generalization 
of standard linear response arguments which gives rise to a ^-diffusion equation including a 
generalization of the Einstein relation for the diffusion coefficient. The equation is similar in 
structure to the porous media equation, but with the important difference that the diffusion 
coefficient depends on the solution of the equation which leads to the fact that, in the absence 
of drift, the diffusion process is classical with mean-squared displacement increasing linearly 
with time. The new g-diffusion equation admits of the well known g-exponential distribution 
which is often postulated as the signature of nonequilibrium systems. 
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